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Abstract 

We estimate the heavy quark effective theory parameter Ai from inclusive semilep- 
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1 Introduction 



As is well known, the heavy quark symmetry breaking parameters Ai and A2 can affect the 
shape of B meson semileptonic decay spectrum substantially. While it is easy to obtain the 
value of A2, the hyperfine splitting term, from mass difference between B and B* mesons, 
it is very difficult to determine the value of parameter Ai, which corresponds to the kinetic 
energy of heavy quark inside a heavy meson. So finding precise value of Ai is very important 
in understanding of heavy meson decay. 

The CLEO collaboration measured the lepton spectrum in the inclusive B Xiu decay 
both by one lepton tagging and by double lepton tagging 0. In single lepton tagging 
data, leptons from secondary charm decay (6 c — > siu) dominate the low lepton energy 
region. These secondary leptons have typically lower energy than the primary ones, because 
they are from c quark decay. To obtain the B —>■ Xiu lepton spectrum in the low region 
from the single lepton tagging data, these secondary leptons must be separated by fitting 
the spectrum with some assumptions and models. 

In Ref. 0, the parameter values of A and Ai were estimated (with fixed value of A2 = 
0.12 GeV^) by using lepton energy distribution of Ei > 1.5 GeV from CLEO data |^ of 
semileptonic decay B Xiu with single lepton tagging. The advantage of using single 
lepton tagging data is small statistical error, though we cannot use low lepton energy part 
of the data {Ee < 1.5 GeV). 

In Ref. ||2| , CLEO collaboration separated B —>■ Xiu from cascade decays of 6 ^ c ^ siu. 
They selected events with tagging leptons of momentum greater than 1.4 GeV, which are 
predominantly from semileptonic decay of one of the two B mesons in an T(45') decay. When 
a tag was found, they searched for an accompanying electron with minimum momentum 
0.6 GeV. The main sources of these electrons are, (a) the secondary lepton from the same 
B, (b) the primary lepton from the other B and (c) the secondary lepton from the other B. 
Lepton from (c) has the same charge as the tag lepton while leptons from (a) and (b) have 
opposite charge to the tag lepton. And leptons from (a) and (6) have different kinematic 
signatures so that their contributions are easy to separate. In the T(45') decay, the B and 
the B are produced nearly at rest. Hence there is little correlation between the directions of 
a tag lepton and an accompanying electron if they are from different B mesons. If they are 
from the same B, there is a tendency for the tagged lepton and the electron to be back-to- 
back. They analyzed the data with double lepton tagging and separated the primary leptons 
from secondary leptons without model dependence. 

In this paper by using the double lepton tagging data, we made a minimum analysis 
to determine value of the parameter Ai. There is one difficulty in fitting for the data, 
as is well known. Since non-perturbative correction up to 1/mf cannot predict the correct 
shape of lepton distribution near the end point, we have to exclude the high energy data 
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points of the distribution. Choosing Eqcd, the maximum lepton energy that one can trust 
the shape of l/rrib expansion, is very important in this fitting. Following Ref . , we choose 
Eqcd = 2.0 GeV. The double tagged data has larger statistical error than the single tagged 
one, but we can use low energy lepton data model-independently. Therefore, this work can 
complement the work of Ref. |^ . 



2 Theoretical details 

Following the heavy quark effective theory (HQET) 0, the mass of a pseudoscalar or a 
vector meson M containing a heavy quark Q can be expanded as 

mM = mQ + A \ , (1) 

where dM = 3, — 1 for pseudoscalar and vector mesons, respectively, and 

Ai = -^{M{v)\K{tDyK\M{v)), (2) 

A2 = TTT^ {M{v)\K^a,,G^''K\M{v)), (3) 

2dMmM 2 

where is the heavy quark field in the HQET with velocity v. Ai parametrizes the mass 
shift due to the kinetic energy of heavy quark inside the meson, and A2 is related to the 
effect of chromomagnetic interaction between heavy quark and light degrees of freedom. In 
case of a. B meson, we can estimate the value of A2 quite accurately from the mass difference 
between B and B* mesons. 

-J- Ax — A2 / j\ 

TUB' = rub + A , (4) 

Irrib 

mB = mb + A , (5) 

zrub 

^(m|. -m|) = A2 + 0(— )^0.12GeV2. (6) 

4 TUb 

Within HQET the lepton spectrum of the semileptonic decays of a 6-fiavored hadron 
{Hb Xqiu) is calculated in the Ref. [^, |^ , and the result is 



and approximately 



dx 



j6 + 5x _{6-4x)e^ (3x - 6)e^ 5{6 - 4x + x'^)e^' 
'\~ (1 ^ + 3(1 

f 5x (2x^-5x)e^ 2(a;3 - Sx^ + 10a;)e'^ 



3 ■ 3(1 -x^^ 
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where 



-L I 







(8) 
(9) 



Kb = -Xi/ml, Gb = 3X2/ml, 

with rriq denoting the mass of the quark q = u, c in the final state, and Vqb is CKM 
matrix element [H. The terms in the second line and third line of Eq. (0) correspond to 
non-perturbative corrections (NP) to leading order Born approximation of the first line. 
Perturbative corrections of the electron spectrum from 6-decay were calculated in various 



references |T0[- The analytic form of order correction is given ^ as 
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dx 



Fi{x,y,e 



with 



where 



Fi{x,y,e'^) = Hi{x,y) + H2{x,y,e^) - H2{x,y,Zm), 



(10) 



(11) 



Hi{x,y) = 2(x - y){xM - X + y)HB{x,y) 
Y ( 

+ —\xi-A + 5x) + yiA -<6x- 5x^) + y'^il + lOx) - 5?/^ 

+ e^[l - 2x + hx^ + t/(5 - 16x) + lly^] + e^(-2 + 6x - ly) + e'^j 

+ lne[a;(-l + 2a;) + ?/(l-4a;) + 2y2 + e2(i + a;-?/)-e^], (12) 



and 



h + zf2 ^ Y,{z){h + zU) ^ Y,{z){h + zU) 



Apz[z) 



8(l-|/)h(^)]Y 8[P3(^)P 
+ \ ln(^)/7 + ^(T^ + [Li2(t/^+(^)) + Li2(«;-(^))]/9 
- yz + Ayp'i{z)Yp{z), 



(13) 



with 



i/B(a:,Z/) = l-ln(l-x)-ln(l-y/x)-2(po/P3>;-l)ln[(l-x)(l-y/x)-e2] 

^Uii, fi-^Uii, fi-^Uii, fi-1^ 

"-^^ LiJl-i^UliJl-^-^^ 



+ ^ 

P3 



Li2(l 
+ Li2 (1 



+ 2rp(F^ + 21ne) + lneln 



p_ / \ xp- 

{w+ -x){w+ -y/x)- 



[x — w^)(y/x — wJ) 



(14) 



And 



= (i_y)3[5a,2^^(5_2x)]-4e^(l-t/)[x^ + ?/(l-5x + 2x^)+?/^(2-a;)] 
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+e^[-a;2 + y{-l + 6x - 3x^) + y^{-3 + 2a;)], (15) 

/2 = _(i_y)[5a;2 + y(5 + i8x + 3x^)+?/^(3-2x)] 

+4e2(l - y)[x^ + y{l - x)] + e^[x^ + y{l - 2x)], (16) 

h = {l-yf[-hx^ +y{-h-9,x + x^) + y^] 

+2e^{l - y)\2x^ + y{2 - 6x + x^) + y'^] + e^[x^ + ?/(l - 4a; + x^) + y\ (17) 
/4 = 5a;2 + y(5 + 28x + 12x2)+?/2(12 + 4a;-a;2)-?/2 

+£^[-4x2 + y(-4 - 4x + 2x^) + 2?/2] _ e^(x2 + y), (18) 
/g = _5 + lOx + ?/(5 + 24x + 8x2)+ 1/2(5 -18x) + 3?/2 

+e2[4 - lOx - 4x2 + ?/(-8 + 18x) - Ay"^] + e^(l - 4x + y), (19) 

/e = 5 + 10x-4x2 + ?/(14+10x)-3^/2-2e2(2 + 3x-2^/)-e^ (20) 

/7 = -5 + 4a; - 4x2 + Qyx - y"^ + [4(1 + x) - lOy] + e^ (21) 

/s = x(l - x) +?/(-l + X + x2) - 2y2a; + y3 _^ e2j^l _ x)(x - y) (22) 

/g = x + i/(l-2x)+y2 + e2(x-y). (23) 

All the parameters and the kinematic variables in the above expressions are listed in Ap- 
pendix. 

After using the above all formulae, the electron distribution in semileptonic decay of B 
meson can be written as 

dTtheory ^ ( ( ^\ ( dT\ 

where is leading order Born approximation, is non-perturbative correction 

using the HQET and the perturbative as correction. We define the CKM-matrix 



independent decay rate 



and then, semileptonic decay rate Tsl can be written as 

rSL = lc\Vcb\^ + lu\Vub\^. (26) 



Since | Vubp <^ | V^bp, we can neglect b —' u decay. Integrating over Ee of Eq. (24), we obtain 

(27) 
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^SL = lc\Vcb\ = To 

where zq and zi are defined as 



zo 11 - '^^4^9{e)] + IzoiGb - Kb) - 2z,Gb 
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zo = 1 - 8e2 + 8e^ - - 24e^ Ine , (28) 
zi = (l-e2)\ (29) 
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and g{e) is a complicated function of e, which can be approximated to 

31\ .0 3 



To obtain the mass ratio e = rric/mt,, we use the relation 

mb-rric = {ms - mo) - ^{Xi + 3^2) ( — (31) 

/ Krric mi)/ 

^ [ruB — mo) ± (~ 1%) ~ rriB — mn = 3.41 GeV. 
We note that if we use instead the other relations, e.g. 

rrib — rric = {rriB* — m/).) (Ai — A2) ( ^ ~ rriB* — m£,* = 3.32 GeV, (32) 

2 KfTic rribJ 

or 

rrib — rric = (m^ — rrTc) ( ^ ^ WEb — rno = 3.35 GeV, (33) 

2 Vrric nib/ 

where TFTb = -{mB + SmB*) and rno = -{itld + Sttid*), 

the values of correction would become as large as ~ +(2 ~ 4%) depending on Ai. 

For b —>■ ciu, Figs, l(a-b) illustrate the dependencies of various corrections on rrib and 
Ai. All figures in Fig. 1 are with \Vcb\ = 0.04. The value of nib determines mainly the 
overall size of decay width, while other parameters determine the shape of the distribution. 
As can be seen from Figure 1(a), we find that the dependence of F^j;^ on rUb is rather weak 
on the contrary to the naive estimation of Tsl oc ml, and is very sensitive to the quark mass 
difference {mb — nic) |T^. Note also that the shapes of Born approximation and perturbative 
correction are almost insensitive to the value of mf,, while non-perturbative correction is 
quite sensitive to both mdmb and Ai. 



3 Results and discussions 



To compare CLEO data with theoretical calculation, we use minimum method with 



X 



E 

Ei<EQCD 



S(E,)-Ft'^--y(E,;e,Ai) 



(34) 



where B{Ei) and cr{Ei) are experimental data of differential branching ratio and error at 
lepton energy Ei, and F^^'^°'^^ {Ef, e, Ai) is theoretical prediction at Ei as a function of param- 
eters e = m^/mb and Ai. We normalized the decay distribution to have branching ratio of 
10.49% as in Ref. ||, B{B Xiu) = (10.49 ± 0.17 ± 0.43)%, i.e., 



0.1049 



SL 



dr 



theory 

dEf 



(35) 



Ee—Ei 
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We note that, because of exact cancellation between Tsl and T theory, F**^'^™^ is independent 
of Fq, and therefore independent of \Vcb\ and m^. i^t^^^^^'y is only indirectly dependent on 
rrib through the definition of x in Eq. (8). Following Ref. |11|, we use the b quark mass, 
TUb = 4.8 ±0.1 GeV, which has been derived from QCD sum rule analysis of the T system 
|T3[| . For as, we use = 0.22, as in Ref. [Q. 

We here comment on determination of Eqcd, which is crucial for this analysis. As we can 
see in the Figs. 1(b) and 2, non-perturbative corrections are significant only in large electron 
energy region {i.e. i?^ > 1.5 GeV). Therefore, we have to include as many data points up 
to Eqcd, in which theory can give correct shape of lepton energy distribution. Otherwise, 
we cannot fully see the effect of non-perturbative correction which determines the value of 
Ai. However, if we include the data points over Ei > Eqcd, the result will be meaningless 
because the shape of the lepton energy spectrum is not reliable above Eqcd region. The 
numerical value of Eqcd can be estimated from the value of 0: For h ^ u decay with 
m„ = 0, Eqcd ~ 0.9-mb/2 ~ 2.15 GeV. For b ^ c decay, smaller smearing range is required 
near end point, and 

Eqcd ~ 0.9 ■ {ml - ml)/2mb ~ 2.0 GeV. (36) 

Since we are dealing only with lepton energies less than 2.15 GeV, we neglect h ^ u decay 
and set 0.6 < Ei < 2.0 GeV. 

We tabulated the results in Table 1. Since we fixed (m^ — m^), changing the value of 
rrtb means changing mass ratio e = rric/mb together, and this mass ratio affects the results. 
The values of A are determined from the mass relation Eq. (Q). All values of Ai in Table 1 
are much larger than the value in Ref. which is Ai = —0.19 ± 0.10 GeV^, or the values 
T3, O] which are ~ —0.1 GeV^, but consistent with JTBl O, IT3| which are in the range 



m 



-0.4 ~ -0.7 GeVl 

The values of Ai show significant dependencies on the input value of rUb, but still each 
value is consistent within la error range. As explained before, this large sensitivity comes 
from mass ratio rric/mb. Indeed, for rrib = 4.8 GeV, if we change the value of {rub — rric) 
to 3.35 GeV, i.e. ml/ ml = 0.091, then Ai = -0.69 ± 0.22 GeVl Changing the value 
of mb with fixed ml/ml = 0.084, we obtain Ai = -0.55 ± 0.18 GeV^ for mb = 4.7 GeV 



Table 1: Results of the fitting with mb = 4.7 ~ 4.9 GeV and the fixed mb — mc = 3.41 GeV. 



rrib 


= ml/ ml 


Ai 


A 


4.7 GeV 


0.075 


-(0.45 ±0.19) GeV^ 


0.57 ±0.018 GeV 


4.8 GeV 


0.084 


-(0.58 ±0.23) GeV^ 


0.46 ± 0.022 GeV 


4.9 GeV 


0.092 


-(0.70 ± 0.27) GeV^ 


0.34 ± 0.026 GeV 
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and Ai = —0.52 ± 0.29 GeV^ for mt, = 4.9 GeV, which are very similar to the case with 
rrif, = 4.8 GeV, as shown in Table 1. 

Once we know the parameter values m^, Ai, we can extract \Vcb\ from the relation 



2 B{B Xiu) 



W = ^ (37) 



For tb, we averaged Tb± and r^o from Particle Data Book |19 



This gives the value 



Tb± = (1.62 ± 0.06) X IQ-'^sea 
tbo = (1.56 ± 0.06) X lO^^^sec 



iKfel = 0.041 ±0.002 , (38) 



where the error includes the errors from semileptonic branching ratio of CLEO data, B meson 
lifetime, uncertainties from Ai and b quark mass. This result is consistent with CLEO result 
with ISGW model which is iKbl = 0.040 ± 0.001 ± 0.002 p[, and with recent Particle Data 



Book resuh = 0.0395 ± 0.0017 |9 



Figure 2 shows the best fit result of differential branching ratio compared with CLEO data 
as a function of charged lepton energy, with rrib = 4.8 GeV and Ai = —0.58 GeV^. It shows 
the relative size and shape of the various corrections for rrib = 4.8 GeV. Non-perturbative 
term is about ~ —4.5% and perturbative term is ~ —12% from leading approximation. 
From these facts, it is clear that non-perturbative correction determines the shape and not 
much effect on total decay rate, while perturbative term has little effect on the shape but its 
contribution on the total decay rate is quite large. Fitting with data points between 1.0 GeV 
<Ei< 2.0 GeV, we obtain Ai = -0.57 ± 0.19 GeV^ and A ~ 0.46 for = 4.8 GeV, which 
are almost the same with the results from 0.6 GeV < Ei < 2.0 GeV. Finally we note 
the dependence on is very weak. Changing to 0.35 we get Ai ~ —0.54 GeV^ and 
A ~ 0.46 GeV for nib = 4.8 GeV, which are almost same values in Table 1. 



We finally note that recently CLEO collaboration measured |2^ the first and the second 
moments of the hadronic mass-squared distribution in the inclusive decay B Xciv and 
also made a preliminary determination of the first and the second moments of the lepton 
energy distribution from the spectrum in Ref. 0. Using those four moments, they ob- 
tained the values of A and Ai, but there appeared to be inconsistencies in the results which 
suggests either experimental error or problems in the HQET. However, if we consider only 
the moments of the lepton energy distribution, the preliminary CLEO analysis pO[ gives 



Ai ~ —0.75 ± 0.20 GeV^, which is rather in a good agreement with our results. 
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A Appendix — Kinematic Variables 

In reference 0, the kinematic variables are defined as: 

• b, q, Gji,^ : four-momenta of the 6-quark, hghter quark, gluon, lepton, neutrino. 

• P = q + G,W = i + i>: four-momentum of the quark-gluon system and the virtual W 

• Ac stands for the scaled gluon mass (Ac = fnc/mb <^ e) 
The scaled masses and lepton energies 



~ rUb \ 
vary in the region 



< X <XM = l-e^ (40) 

< y <ym = x{xm - x)/ (1 - x) (41) 

{e + \Gf< z <zm = {l-x){l-y/x). (42) 

Frequently used kinematic variables which characterize the quark-gluon system are 



1/2 



Pq{z) = \{l-y + z), 
psiz) = ^[l + y^ + z^-2{y + z + yz)] 
P±{z) = Po{z)±P3{z), 



and similarly for the virtual W 

wo{z) = l{l+y-z) 



wsiz) = m + y^ + z^ -2{y + z + yz)] 



1/2 



w±{z) = wo{z) ± W3{z) 

iln^ = ... ^ 

2 w^{z) y/z 



YUz) - iln^ = ln^. (44) 



For G = 0, which implies z = e^, the abbreviations 

^o=Po(e^), P3=P3(e^), etc., 
Wo = wo(e^), W3 = w^3(e^), etc. 

will be useful. Polylogarithms are defined as real functions, and in particular 



(45) 



rx (H 

hUx) = - — ln|l (46) 
Jo t 
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Figure captions 



Fig. 1 Contributions of each terms in lepton spectra of 6 — > civ decay. In all figures, \Vcb\ = 
0.04. (a) Born approximation and perturbative Ug correction with mb = 4.7 GeV (solid 
line), rub = 4.8 GeV (dashed line), = 4.9 GeV (dotted line) and as = 0.22. (b) 
Non-perturbative correction with Ai = 0.3 GeV^ (solid line), Ai = 0.4 GeV^ (dashed 
line), Ai = 0.5 GeV^ (dotted line). A2 and nib are fixed with the values A2 = 0.12 GeV^ 
and TUb = 4.8 GeV. 

Fig. 2 Best fit result for nib = 4.8 GeV with Born approximation (long dashed line), non- 
perturbative correction (short dashed line), perturbative correction (dotted line) and 
sum of the all (solid line). Dots with error bars represents CLEO data. Parameter 
values are Ai = -0.58 GeV^ A2 = 0.12 GeV^ and a, = 0.22. 
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Figure 1: (a) 
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Figure 1: (b) 




13 



